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We report on a theoretical study about the magneto-dipolar coupling and synchronization be¬ 
tween two vortex-based spin-torque nano-oscillators. In this work we study the dependence of the 
coupling efficiency on the relative magnetization parameters of the vortices in the system. For that 
purpose, we combine micromagnetic simulations, Thiele equation approach, and analytical macro¬ 
dipole approximation model to identify the optimized configuration for achieving phase-locking 
between neighboring oscillators. Notably, we compare vortices configurations with parallel (P) po¬ 
larities and with opposite (AP) polarities. We demonstrate that the AP core configuration exhibits 
a coupling strength about three times larger than in the P core configuration. 


I. INTRODUCTION 

In the last decade great attention has been drawn 
to the phase-locking phenomena of spin-torque nano¬ 
oscillators (STNOs)^“^^. STNOs are anticipated to 
be promising devices for sub-micron scale microwave 
synthesizers because of their high emission frequency 
tunability^'^^^®. However, an important issue of such 
devices regarding their practical realization is their low 
output oscillation power and low spectral stability. A 
possible solution to these issues could be the synchro¬ 
nization of a few Synchronization 

between multiple auto-oscillators can also be useful in the 
framework of developing associative memories architec- 
tures^®^^^. Previous studies reported on synchronization 
of STNOs interacting with others via spin waves^®’^®’®^, 
exchange coupling®, electric currents®’^"^’®®, noisy current 
injection^®, or via magneto-dipolar interaction^^’^^’®^“^®. 

Among the various synchronization mechanisms, 
magneto-dipolar coupling is inherent and efficient as em¬ 
phasized in our previous works^^’^^ but also in refs®^“^®. 
In the present study, we focus on the magnetodipolar 
interaction between two vortex based STNOs. 

Single magnetic vortices in cylindrical dots are char¬ 
acterized by two topological parameters®^. Chirality (C) 
determines the curling direction of the in-plane magne¬ 
tization, such that C = -1-1 (resp. C = —1) stands for 
counter-clockwise (resp. clockwise) direction. The orien¬ 
tation of the vortex core magnetization is described by 
its polarity (P), which takes a value of P = -1-1 (resp. 
P = — 1) for core magnetization aligned (resp. anti¬ 


aligned) to the out-of-plane (z) axis. The relative con¬ 
figuration of two interacting vortices can then take four 
non-equivalent states, with identical/opposite chiralities 
and identical/opposite polarities. 

In a previous work^®’^^, we studied the capability of 
two vortex-based STNOs to synchronize through dipolar 
coupling. In this first approach, we have only considered 
the case of two vortices with identical polarities and chi¬ 
ralities, and already demonstrated the possibility to ob¬ 
serve synchronization. In this new study report, we show 
that changing the relative polarity and chirality param¬ 
eters will strongly modify the interaction between the 
auto-oscillators and may strongly modify the efficiency 
of synchronization. We conduct a numerical study in 
which we investigate the synchronization properties for 
selected combinations of vortex parameters, aiming at 
sorting the best combinations of the (C, P) parameters to 
achieve synchronization. We also consider two different 
electrical connections for the current injection i.e, paral¬ 
lel and series connections, corresponding respectively to 
current flowing in the same or opposite direction in the 
two STVOs. 


II. PRESENTATION OF THE SYSTEM 

The studied system consists of two circular nanopil¬ 
lars with identical diameters 2R — 200 nm, separated by 
an interdot distance L (see Fig. 1). Each incorporates 
a Permalloy free magnetic layer (Mg = 800 emu/cm®, 
A = 1.3 X 10“® erg/cm, a = 0.01) with thickness h = 10 


2 


nm, separated by an intermediate layer (non-magnetic 
metal or tunnel barrier) from a polarizing layer with 
perpendicular magnetization. Considering their dimen¬ 
sions, each free layer has a magnetic vortex as its rem¬ 
nant magnetic configuration. The vortices parameters 
will be referred as Pi ,2 and Ci ^2 for 1 st and 2 nd pillar. 
The polarizing layers, whose magnetizations are identi¬ 
cal and oriented along z, will be considered in simulations 
only through the corresponding current spin polarization 
Pzi = Pz 2 = Pz = +0.2. The gyrotropic motion of a 
vortex core can be driven by spin transfer torque action, 
by flowing current above a threshold amplitude through 
each pillar ; In our case, a current density J = 7 x 10® 
A/cm^ (doc = 2.2 mA). Yet, the current sign in each 
pillar has to be chosen so that liPiPz < 0 to ensure 
self-sustained oscillations^^The core polarity of each 
vortex then defines its gyration direction^® (see Fig. 1). 
Indeed, when Pi = +1 (resp. Pi = —1) the vortex core 
circular motion is counterclockwise (resp. clockwise). 


config. 
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Table I. Studied configurations with their respective signs of 
the vortex parameters {Ci, Pi), current density (Ji), and cur¬ 
rent spin polarization (pzi)- 
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Figure 1. (Color online) Schematic illustration of the studied 
system composed of two magnetic dots, each in a magnetic 
vortex conhguration. Vortices cores are shown by small green 
disks. The vortex core positions are given in polar coordi¬ 
nates, i.e. {Xi, (fii) and {X 2 , (^ 2 ), respectively. The blue (resp. 
red) arrows show the core up (resp. down) gyrotropic motion 
sense. 


III. POSSIBLE CONFIGURATIONS 

We then consider six possible configurations for which 
self-sustained oscillations are achieved in both pillars, de¬ 
scribed in Table I. Note that Parallel cores (Pc) config¬ 
urations correspond to vortices moving in the same di¬ 
rection, whereas Anti-Parallel cores (APc) configurations 
correspond to vortices moving in the opposite directions. 

Considering the different configurations displayed in 
Tab. I, the electrical connection must be adapted ac¬ 
cording to the relative vortex core polarities in order to 
fulfill the condition hPiPz < 0 to ensure self-sustained os¬ 
cillations^^As a consequence, the Pc configurations 
must be alimented using the parallel connection to en¬ 
sure same current sign in both pillars (see Fig. 2a). On 
the contrary, the APc configurations have to be supplied 
with a series connection to ensure opposite current signs 
(see Fig. 2b). 


Figure 2. (Color online) a) (resp. b)) Illustration of the DC 
supplied current for the P (resp. AP) configurations showing 
the parallel (resp. series) connection. 


IV. MACRO-DIPOLE ANALYTICAL MODEL 


To get some insights for the origin of the dependence 
in effective coupling with vortices configuration, we con¬ 
centrate in this section on an analytical model based on a 
macro-dipole approximation. The dipolar energy (Wint) 
between two magnetic dipoles /xi and /i 2 is then given by 
the following equation (in CGS units): 


VFint 


( 3 (mi ■ ei2)(M2 ■ 612) - Ml ■ P2) 

IID12P 


( 1 ) 


where D 12 is the vector between the two dipoles and 612 
is a unit vector parallel to D 12 . 

Considering two planar dipoles induced by the off- 
centered vortices in the framework of the two vortex 
ansatz^^ (TVA): ^1 = crCiAi)—sin ((/ji) , cos ((/^i)) and 
^2 = o-C' 2 A 2 (-sin((/? 2 ) ,cos((p 2 )), where a = ^MsV/R, 
^ = 2/3, V = T:P?h . For D 12 = (d, 0), where d = 2R + L 
is the inter-dipole distance along a;-axis, and using equa¬ 
tion ( 1 ) one obtains: 


Wint = -C'iC'2^AiA2(cos((/?i - (P2) - 3 cos((/ 3 i + ^32)) 

( 2 ) 


where pi = PitUi. 

To illustrate the different situations, we consider syn¬ 
chronized oscillations in the two relative polarities config¬ 
urations. For two vortices with same core polarity (Pc), 
gyrating in identical directions at the same frequency 
+1 — +2 ~ 0 and (fi -\- p 2 ~ 2wo: so that equation (2) 
gives: 




int = -C'iC 2 ^AiA 2(1 - 3cos (2a;ot)) (3) 
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In contrast, for vortices with opposite polarities (APc), 
gyrating in opposite directions, Lpi-\-Lp 2 « 0 and Lpi—Lp 2 « 
2 a;o, so that one obtains: 

= -CiC 2 ^XiX 2 (cos (2a;ot) - 3) (4) 

Equations (3) and (4) show that for a given vortex 
gyration frequency loq the coupling energy W-mt oscillates 
at twice the frequency (2wo). In the Pc case (see blue 
curve in Fig. 3) it oscillates with a large amplitude and 
a small mean value, whereas in the APc case (see red 
curve in Fig. 3) it oscillates with a large amplitude and 
a smaller mean value. 



Figure 3. (Color online) a) Dipolar energy (Wlnt) evolution 
of two interacting vortices modelled as macro-dipoles and os¬ 
cillating at the same frequency. The blue curve corresponds 
to the identical polarities (Pc) case and the red curve cor¬ 
responds to the opposite polarities (APc) case. The dashed 
colored lines represent the corresponding mean value of the 
coupling energies (Wint)- 


When not synchronized, the two vortices will feel two 
oscillating components of the magneto-dipolar interac¬ 
tion i.e., one at low frequency and one at high frequency. 
The latter one will average out and have negligible in¬ 
fluence on the phase locking features, while the low fre¬ 
quency term will be responsible for the synchronization 
phenomenon. The effective coupling coefficient /Teff can 
be identified writing (VFint) = HeffCiC 2 XiX 2 for the 
mean coupling energy and gives the following results for 
the Pc and APc relative vortex core polarity configura¬ 
tions: 


Meff = - 


= 3 


2 <73 


2 # 


Synchronized states correspond to a minimization of 
the average interaction energy. As illustrated here, rela¬ 
tive polarities and chiralities signs influences the sign of 


Wint- As a consequence, these relative parameters also 
define the phase relationship achieved when synchroniza¬ 
tion occurs. The later considerations are illustrated in 
figures 4 and 5. 

From this study, we then conclude that the effective 
coupling coefficient is predicted to be three times stronger 
when polarities are opposite (APc) than when polarities 
are identical (Pc). Concurrently, the high frequency os¬ 
cillation of interaction energy is three times larger in Pc 
polarities configuration as compared to APc case. While 
this indicates that APc is the optimal configuration for 
synchronization, we must note that this second contribu¬ 
tion may affect the locking phenomenon. 
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C, = -l =G 
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Figure 4. (Color online) Schematic illustration of the coupled 
dynamics for the Pc configurations, i.e. where Pi = P 2 - 


V. THIELE ANALYTICAL APPROACH 


The spin transfer induced gyrotropic vortex dynamics 
can be described by the Thiele equation^^’^^’"^®’^®: 

P,G, X X, -h D, • X, - UXi, Ci,Ji)X, 

(5) 

-FpT(p,„ PO - Fi„t(X„ P,P„ C,C,) = 0 

where = —Gz is the gyrovector with G = 27rMs/i/7 

and Di = arjiG is the damping coefficient with iji = 
O.Sln (i?i/(2Lex)) + 3/8,^^ where Lex is the exchange 
length. For each pillar, the vortex frequency^^’^® is given 
by the ratio between the confinement coefficient ki and 
the gyrovector Wop = ki/G with: 


fc,(x„ Q, j,) = kr+krc^J^+ik'r+k'^^c,j,) 
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Figure 5. (Color online) Schematic illustration of the coupled 
dynamics for the APc configurations, i.e. where P\ = —P 2 - 


where A™® and (resp. and k'f’^) corresponds 
to the magnetostatic (resp. Oersted field) contribution. 
The Oersted contribution will increase the vortex core 
gyration frequency if the vortex chirality is along the 
same direction as the Oersted field {CiJi > 0), and re¬ 
spectively decrease the frequency otherwise^®. Gyration 
amplitudes will also be affected by such interaction with 
Oersted field. To maximize the symmetry of the system 
and avoid that the Oersted contribution brings an offset 
between the two STNOs’ frequencies, we find that the 
condition C 1 J 1 C 2 J 2 > 0 should be ensured (correspond¬ 
ing to identical Oersted contributions in both pillars). 
This excludes configurations Pc3 and APc3 from Tab. I 
from being optimal configurations for synchronization. 

The fourth term in Eq. (5) is the spin transfer force, 
which for the case of perpendicularly uniform magnetized 
polarizer writes^^: 

rfT = K(X,xz) (7) 

where k = njajMsh is the effective spin torque effi¬ 
ciency on the vortex and aj = lipzPJ/{2\e\hMs). In this 
study, we chose to neglect the field-like torque (FLT) 
contribution. While the FLT is negligible in case of a 
metallic intermediate layer, its amplitude can however 
reach a significant fraction of the Slonczewski torque in 
case of a magnetic tunnel junction. However, micromag- 
netic simulations computed with a FLT contribution of 
10% (typical) of the magnitude of the Slonczewski term 
showed no significant influence on the gyrotropic dynam¬ 
ics. A last term accounts for the interaction dipolar force 
between the two neighbored vortices: Fintji(Xi 2 ) = 
-d{Wint) (Xi,X2)/5Xi, 2 = -ClC 2 MeffX 2 ,l,’ where 
is either or depending on P 1 .P 2 sign. 

The system of coupled equations for the vortices core 
motion given in Eq. (5) provides a dynamical descrip¬ 


tion of the phase locking between the two cores. We 
introduce the two variables 'F = Pipi — P 2 P 2 and e = 
{Xi — X 2 )/{Xi + X 2 ). Following the methodology de¬ 
scribed it Belanovsky et al. , by linearizing the system 
around equilibrium trajectories, we obtain a linear set of 
equations describing the evolution in time of the relative 
phases and amplitudes: 

I e = -2ar] -F e - ^4- (8a) 

W =-4 £-|-2q!77^4' (8 b) 

where rg = Xq/R is the normalized average gyration 
radius and a = k'^^/kms = 1/4. The two equations (8a) 
and (8b) are linear and their eigenvalues are 

Al 2 = -ar]UJoarl±^a'^ifujl^arlY + 4^ - 

In the case of periodic solutions, the phase-locking dy¬ 
namics is characterized by a phase-locking time (r) and 
a beating frequency (11) that can be written as: 

{ 1/t = —arjuJoarQ (9a) 

fl^ =-(aTywoaro)^ - 4 - d^woorp (9b) 

In next section, we propose to realize micromagnetic 
simulations®*^, from which H and r will be extracted from 
the phase-locking dynamics. The effective coupling coef¬ 
ficient in each configuration peS will then be derived for 
each considered configuration by simply reverting equa¬ 
tions (9a) and (9b): 

PesiT,n) = ^ (l/(ra? 7 ) - \/l/(ra;? 7)2 - fl(L)2^ . (10) 

These micromagnetic simulations represent a more re¬ 
alistic picture of the coupled system as it takes into ac¬ 
count the non-punctual geometry of the magnets as well 
as the full current induced Oersted field contribution, in¬ 
cluding cross-talk between nano-pillars. 


VI. MICROMAGNETIC SIMULATIONS 

We first compare the results of micromagnetic simula¬ 
tions obtained for the two cases Pci and APcl with a sep¬ 
arating distance between nano-pillars L = 50 nm. The 
evolution of radii and dephasing parameter 'F is shown 
in figures 6(a) and 6(b) respectively and some numerical 
values are given in Tab. 11. These results first confirm 
that phase-locking is achieved in both configurations. For 
both configurations, self-sustained unlocked oscillations 
in each pillar start at the same frequency but start with 
a random phase shift, and then converge towards phase- 
locked regime in very close phase-locking times (r). In 
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01 C, = -1 P„sym1 C 2 = -1 gj] C, = -1 AP_sym1 C 2 = +1 



P, = -1 P2 = -1 Pl = -1 P2 = +1 




Figure 6. (Color online) a) (resp. b)) Vortex cores orbital 
radii and phase difference ^ = ipi — ip 2 (resp. sum = 
‘fii + ipa) obtained by micromagnetic simulations for the Pci 
(resp. APcl) configuration where L = 50 nm. As seen in 
top schematic illustrations the chirality (C) is opposite to the 
Oersted field (OH) in all the dots (C/OH[ti]). 


config. 

/ (MHz) 

Xoi (nm) 

X 02 (nm) 

vk 

Pci 

468.80 

63.59 

63.59 

^ 0 

APcl 

470.57 

64.46 

64.46 

^ 0 

APc2 

497.75 

44.51 

44.51 

0 

APc3 

476.31 

65.59 

40.59 

—^ TT 


Table II. Numerical values of parameters extracted from mi¬ 
cromagnetic simulations: / is the common oscillation fre¬ 
quency, Xoi (resp. A 02 ) is the left (resp. right) dot vortex 
steady-state radius, and ^ is the dephasing parameter. 


config. 

T 

Q 

fJ-es/G 

(IFint) 

(ns) 

(MHz) 

(MHz) 

(xlO-^'^ erg) 

Pci 

82.78 

40.136 

19.7 

-22.75 [-27.08] 

APcl 

71.20 

67.380 

49.2 

-58.31 [-64.23] 


Table III. Numerical values of r, fl, fign/G, and (Hint) ob¬ 
tained after combining micromagnetic simulations and our 
Thiele equation approach for Pci and APcl configurations. 
The last column containing the mean interaction energy com¬ 
puted by Eq. (10) shows also the numerical evaluation using 
Eq. (11) inside brackets. 


their phase-locked state, both vortex cores oscillate with 
identical radii. 

The phase dynamics obtained by micromagnetic sim¬ 
ulations are fitted to fk = sin(nt -|- ipo) to extract 

n the beating frequency and t the convergence time for 
phase-locking (see Tab. III). The effective coupling val¬ 
ues for L = 50 nm are then deduced: fies/G = 19.7 
MHz for Pc configuration, and /ieff/G = 49.2 MHz for 
the APc one. The coupling strength then appears to be 
stronger in the AP configuration (~ 2.5 x) as expected 
from macro-dipole model. 


The results for the ”APc2” and the ”APc3” configu¬ 
rations for L = 50 nm are shown in Fig. 7. Again in 
both cases phase-locking is achieved. In the symmetric 
case APc2, for which both chiralities are parallel to the 
Oersted field, starting frequencies are again identical in 
each pillar, whereas it is not the case for APc3 config¬ 
uration, in which symmetry is broken by the Oersted- 
field being opposed to chirality in only one pillar. In 
the latter case, the two auto-oscillators have to adapt 
their frequencies to achieve synchronization to a com¬ 
mon frequency /i = /2 = 476.31 MHz, by shifting their 
amplitudes accordingly. As highlighted previously, the 
micromagnetic simulations confirm that the equilibrium 
phase shift changes from jlkl = 0 to |'I'| = tt when the 
sign of respective chiralities sign(GiG 2 ) changes. 


HI C, = +1 AP_sym2 C 2 = -1 [3] C, = -1 AP^asym C 2 = -1 





Figure 7. (Color online) Vortex cores orbital radii and phase 
sum T = ifii + (fi 2 where L = 50 nm for configurations: a) 
”APc2” where both the vortex chiralities are aligned with 
the Oersted field (C/OH[tt]), b) ”APc3” where the vortex in 
the right (resp. left) dot has an anti-aligned (aligned) chi¬ 
rality with the current induced Oersted field C/OH[74] (resp. 

(C/OH[tt])). 


VII. NUMERICAL APPROACH 

To investigate further the difference in coupling 
strength between Pc and APc configurations, and vali¬ 
date the macro-dipoles approach, a more precise numer¬ 
ical calculation of the dipolar energy is proposed. The 
dipolar interaction energy is here summed up over the 
full magnetization distributions obtained by micromag¬ 
netic simulations. It consists in taking into account all 
the spin to spin i.e., cell to cell, interactions between the 
left pillar and the right pillar as follows: 

Ni N2 

^nr = EE^-bb> (11) 

i=i i=i 

where IFi„t,y = -(3(^i • ■ Sij) - fj,i ■ fij)/\\T>ij\\^. 

Ni (resp. N 2 ) is the number of cells in the left (resp. 
right) dot. 
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As illustrated in Fig. 8 each dot can be seen as com¬ 
posed by two distinct regions. The outer part (OP) and 
the inner part (IP) with respect to the vortex gyrotropic 
trajectory. The OP is a quasi-static region and the in¬ 
ner part can be considered as an oscillating dipole. For 
means of comparison with our analytical model the OP 
region is first neglected. As seen in Fig. 10 the val¬ 
ues of (W;™™) are close to the values of (Wint) obtained 
through the macro-dipoles and Thiele equation approach 
when the OP region is neglected. 



Figure 8. (Color online) Illustration of the in-plane magneti¬ 
zation of the two oscillating vortices. The gray zone represent 
the quasi-static magnetization (outer part - OP, r > Aq) and 
the green zone represents the oscillating magnetization (inner 
part - IP, r < Xq). The dashed lines give the vortex orbital 
movement delimitation where r = Xq. 

Figure 9 shows the results for an edge to edge dis¬ 
tance between two STVOs of L = 50 nm considering the 
IP region only in which we compile the data for APcl 
configuration (red triangles) and Pci configuration (blue 
dots). The dashed lines give the mean value of the in¬ 
teracting dipolar energy . In both cases and as 

expected, the energy Wint™ oscillates at a frequency that 
corresponds to twice the gyrotropic frequency (see Tab. 
II). 

We reproduced the process for several other distances 
between the dots {L = 100,200,500 nm). The evolution 
of the average interaction energy versus L extracted from 
micromagnetic simulations is shown in Fig. 10. For both 
Pc and APc configurations. The agreement between nu¬ 
merical and Thiele-based estimation of interacting energy 
is fairly good, notably for large interdot distance. 

As discussed theoretically in section IV, the macro¬ 
dipole model gives a ratio of 3 between the interaction 
energy for the P and AP core configurations. In contrast 
to that prediction, it should be noticed that for small 
values of the edge to edge interdot distance L the ratio 
between the calculated energy gets smaller (~2.6 for L = 
50 nm). 

The thermal fluctuations were not involved in our 
simulations. Yet, the obtained mean coupling energy 
I (Wint) I can be easily compared against thermal energy 
{k^T). The condition for stable synchronization i.e., 

I (Wint) I > fceT, are then found to be D12 < 400 nm 
in the Pci configuration and D 12 < 550 nm in the APcl 
configuration. Finally, the data is fitted with a decay law 
for Wint as function of D 12 , i.e. D ^2 where a has values 
between 3 and 4 as shown in of figures 10 and 11. 


Num. calc. (P_syml) i -A- Num. calc. (AP_syml) 



For /. = 50 nm 


0 12 3 4 

Time (ns) 

Figure 9. (Color online) Numerical computation of the dipo¬ 
lar energy (W;"™) evolution of two interaction and synchro¬ 
nized vortices extracted from micromagnetic simulations for 
L = 50 nm considering the IP region only (see Fig. 8). The 
blue curve shows the evolution of the in phase oscillating vor¬ 
tices (with parallel polarities, P conhguration) and the red one 
the anti-phase case (vortices with anti-parallel polarities, AP 
conhguration). The dashed colored lines represent the corre¬ 
sponding mean values of the coupling energies ((Win™)). 



D ,2 = 2R + L (nm) 


Figure 10. (Color online) Inter-dot distance {D 12 = 2R + L) 
dependence of the absolute value of the mean coupling en¬ 
ergy (Wint): The case for parallel (resp. anti-parallel) po¬ 
larities using a macro-dipoles and Thiele equation approach 
model [hlled squares (resp. unhlled squares)] and by numeri¬ 
cal dipole-dipole computation of the inner parts (IP, see Fig. 
8) of the vortex cores trajectories [hlled circles (resp. unhlled 
circles)] are shown in blue (resp. red). 


VIII. CONCLUSION 

In conclusion, we performed a comparative study of 
vortices parameters configuration for the synchronization 
of two dipolarly coupled spin transfer vortex-based oscil- 
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By studying different contributions to the coupled vor¬ 
tices dynamics, we have also shown that this configu¬ 
ration matches with feasible experimental configuration. 
Optimal configuration then corresponds to nano-pillars 
connected in series. 

Comparing the computed Wint (IP only) with the ther¬ 
mal energy k^T, one obtains that synchronization can 
be presumably achieved D 12 ^ 400 nm in case of parallel 
polarities configuration (Pc), while £>12 ^ 550 would be 
sufficient in case of anti-parellel polarities. 

As far as phase-locking stability is concerned, we high¬ 
lighted that the dipolar interaction keeps involving strong 
oscillations in the coupling energy even after achieving 
synchronization. These interactions will play against 
synchronization and should decrease the minimum inter¬ 
pillar distance to achieve synchronization. 


Figure 11. (Color online) Inter-dot distance (D 12 = 27? -|- 
L) dependence of the absolute value of the coupling energy 
(Wint): The case for parallel (Pc) [resp. anti-parallel (APc)] 
polarities obtained by numerical dipole-dipole computation 
for the whole volume (IP -|- OP) [filled triangles (resp. unfilled 
triangles)] are shown in blue (resp. red). 

lators. As the major result of this numerical and ana¬ 
lytical study, we demonstrate that the effective coupling 
of two vortices with opposite core polarities and hence 
gyrating in opposite directions is larger than the case of 
identical polarities. 
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